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Abstract In complex transition metal oxide heterostructures of physically dissimi- 
lar perovskite compounds, interface phenomena can lead to novel physical proper- 
ties not observed in either of their constituents. This remarkable feature opens new 
prospects for technological applications in oxide electronic devices based on nm- 
thin oxide films. Here we report on a significant electronic charge and orbital re- 
construction at interfaces between YBa2Cu3C>6 and SrTiOs studied using local spin 
density approximation (LSDA) with intra-atomic Coulomb repulsion (LSDA+C/). 
We show that the interface polarity results in the metallicity of cuprate-titanate su- 
perlattices with the hole carriers concentrated predominantly in the CuC<2 and BaO 
layers and in the first interface T1O2 and SrO planes. We also find that the interface 
structural relaxation causes a strong change of orbital occupation of Cu 3d orbitals 
in the CuC<2 layers. The concomitant change of Cu valency from +2 to +3 is re- 
lated to the partial occupation of the Cu 3fi? 3 .2_ r 2 orbitals at the interface with SrO 
planes terminating SrTiC>3. Interface-induced predoping and orbital reconstruction 
in CuC»2 layers are key mechanisms which control the superconducting properties 
of field-effect devices developed on the basis of cuprate-titanate heterostructures. 
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1 Introduction 

It is well known that the rich properties of transition metal oxides like ferroelec- 
ticity, magnetism or superconductivity are closely associated with the physics of 
c/-orbitals (TJ. In heterostructures of physically different transition metal oxides, the 
local structural deformations and uncompensated charge at the interfaces result in 
properties well beyond a simple combination of the characteristics of their con- 
stituents J2] [3] |4). A prominent example is the titanate superlattice of insulating 
perovskites SrTi03 and LaTi03 where the metallic conductivity is caused by the 
mixed valence (+3/ +4) of Ti and by electronic charge reconstruction at polar in- 
terfaces (5). The new electron states appearing in such heterostructures can be tuned 
by external electromagnetic fields which opens new prospects for engineered oxide 
electronic devices based on nm-thin transition metal oxide films. 

In heterostructures with high-r t cuprate films, the external electrostatic fields 
tune the surface/interface charge density and in this way lead to transitions between 
conducting (or superconducting) and insulating states. Such heterostructures consist 
of YBa2Cu307_g-films grown on SrTiC>3 layers and are of central importance in su- 
perconducting field-effect devices |6] [7] ID [TO) . In the existing literature, the shift of 
T c achieved with increasing electrostatic field is often understood in terms of an 
electrostatic doping concept without detailed consideration of the microstructure of 
the interfaces between YBaiCu30 7 _,5-fiIm and SrTi03 ||7][9]. Several experimental 
facts, however, indicate an interface-related modification of the electronic states in 
the cuprate/perovskite oxide heterostructures. For example, recent studies on under- 
doped cuprate films produced a r t -shift of about 5-15 K, whereas in the overdoped 
films no shifts were observed, a fact, which cannot be explained satisfactorily by 
field-doping [9|. Despite the important role of interfaces in the physical properties 
of such oxide heterostructures, no direct experimental studies of the electronic prop- 
erties of interfaces between copper and titanium oxides are currently available. This 
is partially related to extreme experimental difficulties in distinguishing the direct 
interface- and bulk-contributions in the experimental data J8). On the other hand, 
due to the high complexity of the cuprate and titanate structures, extensive theoreti- 
cal studies of the cuprate-titanate interfaces HT] [T2] [T3] H3] US] E6) are a challenging 
task which requires extremely demanding computational resources lfT7l[T8l . 

The goal of the current project is the development of interface models for mi- 
crostructures of strongly correlated cuprates and titanates. The analyses of novel 
electronic states are based on the density functional theory within LSDA and 
LSDA+C/ implementations. Due to the insufficient structural information about 
atomic arrangements at cuprate-titanate interfaces, the theoretical investigations in- 
clude the theoretical modeling of heterostructures with different structurally com- 
patible surface terminations of SrTiC>3 and YBa2Cu3C>6 and relaxation of interface 
atomic structures. In this context, the main amount of the computer time is used 
for high-performance theoretical engineering of multilayer oxide supercells and for 
structural optimization. 

In our extensive analyses, we find a metallic state in cuprate-titanate interfaces 
which is characterized by significant hole predoping of cuprate films ffTTl . In the 
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field-effect experiments, such an interface-caused predoping occurs even before the 
electrostatic injection of charge and is of essential importance for functional proper- 
ties of superconducting field-effect devices. Furthermore, out results not only allow 
to obtain a clear answer to the question about the electronic charging of the cuprate- 
titanate heterostructures, but also to make important predictions how to increase the 
performance of field-effect devices. Our studies directly show that despite the ex- 
treme computational resources required for optimization procedures, the interface 
structural relaxation is of crucial importance for understanding the physics of tran- 
sition metal oxide heterostructures. In the copper oxide layers of cuprate-titanate su- 
perlattices, we find that the interface-caused local structural deformations strongly 
affect electronic occupancies of c/-orbitals of Cu and can lead to a change of its 
electronic valence state. Due to the central role of c/-orbitals in the magnetic and su- 
perconducting properties of transition metal oxides, the obtained interface-induced 
orbital reconstruction should have dramatic consequences for the functionality of 
field-effect devices based on cuprate-titanate superlattices. 



2 Description of the methods 

Our extensive studies of oxide surfaces are based on the density functional theory 
within the local spin density approximation (LSDA) with intra-atomic Coulomb re- 
pulsion (LSDA+J/) |[T9l . The calculations of the electronic densities of states, struc- 
tural relaxation in terms of minimization of total energy and forces have been per- 
formed using the linearized augmented plane wave method (LAPW) implemented 
in the WIEN2k package (20). The most numerically intensive parts of the package 
code are executed in parallel, namely full potential calculations, iterative solution 
of eigenvalues and eigenvector problems of the generated Hamiltonian matrix and 
calculations of electronic charge density expansions. Parallelization is achieved on 
the k-point level by distributing subsets of the k-mesh to different processors and is 
based on parallelization libraries including MPI, Scalapack, and BLACS. 

The small amount of the computational resources was spent for the development 
of interface-related extensions to the WIEN2K-code which are aimed at specific 
studies of the local electronic structure and the charge density profiles across the 
interfaces in heterostructures. For instance, to obtain the profiles of the hole density 
across the polar interfaces, we had to generate the charge densities in the energy 
window between the Fermi level and the top of the valence bond and to perform a 
subsequent planar integration of these profiles. In this way, not only the charge dis- 
tribution, but also the effective thickness of the interface, where the charge density 
deviates from its usual bulk-type behavior, can be estimated. 

The main part of the resources was absorbed by computationally demanding 
structural optimization and LDA+C/-calculations of the optimized superlattices. 
Specifically, for superlattices of YBa2Cu306/SrTi03 and Cu02/SrTi03 containing 
from 18 to 22 atoms in a supercell, a single run of LDA+U included from 50 to 
80 iteration steps until the final convergence with respect to the electronic charge 
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and total energy could be reached. With 15-20 processors used for the parallel 
^-points calculations, each converged LDA+C/ run usually required up to 25,000 
CPU hours. In addition, in each superlattice, the relaxation of the interface distances 
could be achieved by performing up to 10 single LSDA runs which required about 
10,000 CPU hours. Due to a wide range of superlattices analyzed in our project 
(YBCO/STO-sandwiches, Cu0 2 /STO-slubs, Sr 2 Cu0 2 Cl 2 /STO superlattices), such 
computationally demanding calculations have used about 150,000 CPU hours in the 
first 6 months after beginning the project. 

The further stages of calculations of interfaces require a full optimization of all 
local atomic coordinates in the superlattices and the increase of the number of el- 
ementary unit cells of the constituent compounds in each supercell. Consequently, 
the computational time required for a structural optimization drastically increased 
up to 25,000 CPU hours for a single superlattice, whereas each LDA+U single run 
requires about 35^40,000 CPU hours for a supercell with 33 different atoms. With- 
out the exceptional computational resources provided by HLRB, the performance of 
such demanding calculations would not be a feasible task and could not have been 
achieved in a realistic time period. 



3 Electronic properties of cuprate-titanate superlattices 

The recent X-ray studies of the interface arrangement ETl [8) give clear indica- 
tions of incompletely grown unit cells of the cuprate film on the SrTiOs substrate. 
Whereas the bulk unit cell of YBa 2 Cu30 7 _5 is terminated by CuO chains, respon- 
sible for the charge doping of cuprates, the CuO chains near the interface are miss- 
ing which makes the corresponding unit cell incomplete. This structural rearrange- 
ment should strongly affect the superconducting properties of cuprate thin films. 
Therefore, we focus in our studies on the consequence of the effective 'substitu- 
tion' of electrostatically charged CuO-chains with electrostatically neutral (001) 
layers of SrTiOs. As the superlattices, which result from such a structural rear- 
rangement, will be polar, the overall charge neutrality requires extra hole carriers 
to be distributed near the cuprate-titanate interface. Our results clearly show that 
in YBa 2 Cu30g-SrTi03 heterostructures such interface induced hole predoping may 
drastically change the interfacial electronic properties. 

In YBa 2 Cu306-SrTi03 heterostructures, the chemical bonding at the (001) in- 
terface with the cuprate film will be determined by the first termination plane of 
SrTi03, which can be either Ti0 2 or SrO. Below we demonstrate the consequences 
of interface electronic reconstruction for these two basic cases of SrTi03 termina- 
tion. 
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3.1 Case study: TiO 2 termination at interface 

In the case when the substrate of SrTiC>3 is terminated by a TiO?-layer, the deter- 
mined interface arrangement is typically a stack of . . ./TiC^/BaO/CuC^/Y/CuCh/. . . 
layers. Such structural stacks suggest an interface chemical bonding Ti-O-Cu with 
the oxygen of the BaO-layers shared between the CUO2 and Ti02-planes. From the 
electrostatical point of view, the initial 'bulk-type' electronic charging of the con- 
stituent layers indicated in the left panel of Fig.Q~]would result in 1 extra hole which 
is needed in order to compensate the polarity. In the sandwich-type supercell shown 
in Fig.Q] this compensation leads to a doping of upper and lower symmetric blocks 
. . ./SrO/Ti0 2 /BaO/Cu0 2 by 0.5 holes. 




SrO 



TiO, 

Fig. 1 Scheme of a YBaCuO/SrTiOj-sandwich where the polar interfaces appear due to the in- 
complete Ba^Cu^Og-unit cell with an interface structural configuration shown in detail on the 
right. 

For the YBaCuO/SrTiO sandwich-type supercells, we performed calculations of 
the electronic structure using the SIC variant of the LSDA+f/ method on a 9 x 9 x 1 
k-point grid with U = 8 eV and J = 0.8 eV on the Cu 3d orbitals. The lattice 
constants of the supercell were fixed to the lattice constant a = b = 3.898 A of 
cubic SrTiC<3. The interface distance A between the nearest BaO and TiC»2 planes 
was optimized by minimization of the total energy which leads to the value A = 
1 .85 A between the apical oxygens of BaO and the T1O2 plane. 

From the calculated density of states, shown in Fig. 12 we can identify the metal- 
lic state with hole carriers originating from the oxygen p states. Like in the bulk 
YBa2Cu3C<6, the Cu d x 2_ v i states are separated by a gap of 1.34 eV from the oxy- 
gen 2p states. It should also be noted that for this type of interface arrangement, the 
orbitals Cu d^i_ r i and t2 S are occupied and remain below the Fermi level. Fig. [2] 
demonstrates that the significant amount of hole charge is located in the interfacial 
Ti02, BaO and more distant SrO layers. Also, a substantial amount of hole charge 
is distributed over the CUO2 planes. 
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Fig. 2 Density of states of the SrTi03/YBa2Cu20(,/SrTi03-sandwich calculated within the 
LSDA+C approach with U = 8 eV and J = 0.8 eV for the electrons in Cu 3d orbitals. The zero of 
energy is at the Fermi level. 



More detailed information about the distribution of the charge-compensating 
hole density is presented in Fig. [3] Through an integration of the hole carrier den- 
sity along the interfacial z-direction we obtain that the major part of hole carriers 
is distributed along the BaO (25%) and the first Ti0 2 (48%) and SrO (12%) layers. 
Furthermore, about 5% of the total interface hole charge 0.5 e is located in CuOo 
planes. This suggests a finite metallic conductivity in the titanate, BaO, and copper 
oxide planes. It is noteworthy that the structural relaxation for this type of interface 
leads to a reduction of the hole density in the BaO plane and to its slight redistri- 
bution within the interface SrO and more distance planes, whereas the hole density 
and orbital occupation in the CUO2 planes remains almost unaffected. 
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Fig. 3 Hole density distribution near the YBa2Cu20(,/SrTi03-interface for optimized (A = 1 .85 A) 
and unrelaxed (A = 1.94 A) cases. The position z = is at the bottom Ti02-plane of the 
SrTi0 3 /YBa 2 Cu 2 6 /SrTi0 3 -sandwich. 
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Fig. 4 Scheme of a polar CuC^/SrTiOj-superlattice where a STO-layer is terminated by a SrO- 
plane. The right side shows a structural configuration which appears at the interface. 



Our results for the hole charge in the copper oxide planes clearly show that, apart 
from chemical doping, the interface polarity is another decisive factor that modu- 
lates the doping level in the cuprate films. In superconducting field-effect devices, 
operated by electrostatic charging, such initial interface-induced predoping levels 
may have dramatic consequences for their performance. From the point of view of 
superconductivity, the most important feature is the shift of T c with electrostatic 
doping which will be directly affected by the obtained predoping (x = 0.025) in 
the copper oxide film. Moreover, in other structural configurations at YBCO/SrTiO 
interfaces, much higher hole predoping levels and dramatic changes of electronic 
occupancies of Cu 3d orbitals in CUO2 planes can appear as a result of electronic 
and orbital interface reconstruction. 



3.2 Case study: SrO termination at interface 

To analyze the structurally different SrO termination of the substrate SrTi03, we 
consider a model case in which a copper oxide plane is directly deposited on SrTiOs 
(Fig. 3J. The direct deposition of the single Cu 2+ 4 ~ plane on the electrostati- 
cally neutral titanate layer would require two extra holes to maintain the overall 
charge neutrality. To achieve such an extremely high doping level, interface elec- 
tronic reconstruction is inevitably required. Apart from the electronic mechanism, 
other forms of interface reconstruction like oxygen vacancies or cation intermixing 
could modify the chemical composition. However, it is still instructive to enforce 
atomically flat and stoichiometric surfaces in order to study comprehensively differ- 
ent mechanisms of the electronic reconstruction. 
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In our work, in order to focus on the effect of electronic reconstruction, we have 
introduced a decoupling vacuum layer of 13 A thickness between the CuC>2 surfaces 
in the slub geometry shown in Fig. [4] We have also performed an optimization of 
the superlattice structure where a relaxed interface distance A = 1 .83 A between the 
CuC>2 and SrO corresponds to a minimum of the total energy. 



4 

< 3 
2 

JZ 

c 1 





Cu0 2 




I 1 I 1 

A=2.2 A 






— A=1.83A ; 




SrO 


Ti0 2 SrO 



2 4 6 8 

z [001], A 

Fig. 5 Hole density distribution in the interface planes of the Cu02/SrTiC>3 superlattice. Here z = 
corresponds to the location of the lower Cu02-plane. 



In contrast to the case of T1O2 termination, the effect of the structural relaxation 
on the distribution of the interface-predopedhole density is crucial for the CuCh/SrO 
interface (Pig.[5j. To demonstrate this fact, we show in Fig.[5]the distribution of hole 
carrier density for the optimized (A = 1 .83 A) and unrelaxed (A = 2.2 A) interfaces. 
The density plots across the interface are obtained by the integration of generated 
charge density in the energy window between the Fermi level and the top of the O 
2p bands. As can be seen in Fig. [5] the relaxation of the interface distance A leads to 
a dramatic decrease of the hole density «/, in the Q1O2 layers. The integration of «/, 
along the interfacial z direction shows that the total predoped hole charge Nh accu- 
mulated predominantly in the 2p orbitals of O and hybridized with 3d Cu orbitals, 
amounts exactly to 2 holes in the unrelaxed (A = 2.2 A) supercell. In contrast to this, 
in the relaxed structure the density integration gives Nf, = 1 hole for the total hole 
charge. Such a strong difference between N], in relaxed and unrelaxed cases appears 
due to the additional orbital mechanism of the polarity compensation. It becomes 
important due to the interface structural relaxation. 

The new mechanism of the orbital reconstruction becomes evident when we com- 
pare the partial densities of states of 3d orbitals calculated for optimized and unre- 
laxed supercells demonstrated in Fig. [6] The figure shows that the 3d x 2_ji and 3d xy 
orbitals of Cu are affected only slightly by the decrease of the [Cu02-SrO] distance 
A . Specifically, the orbitals 3d xy are fully occupied and located well below the Fermi 
level independently of A . Furthermore, a significant amount of holes is found in the 
states 3d x i_ y i which are split due to the charge transfer gap of about 1.5 eV. The 
charge transfer gap becomes narrower with the decrease of A . The most striking ef- 
fect, which appears for the relaxed 4 = 1.83 A, is the reconstruction of Cu 3c/ 3 _2_ r 2 
orbitals. While for A = 2.2 A both Cu 3c/ 3 _2_ r 2 orbitals are located about 4 — 5 eV 
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Fig. 6 Projected Cu 3d orbital density of states at the interface between CuOi and SrTiOj termi- 
nated by SrO (LSDA+[/ studies). The zero of energy is at the Fermi level. The left and right panels 
correspond to the case of unrelaxed (A = 2.2 A) and optimized (A = 1.83 A) interfacial distances. 



below the Fermi level, the decrease of A leads to their splitting. This splitting re- 
sults in the "displacement" of one d 3 2_ r 2 orbital (spin down case in Fig.|6]l to about 
1 eV above the Fermi level which implies that this orbital becomes empty. Such a 
strong change of the occupation of d^ z i_ r i corresponds to an increase of the valency 
of Cu from +2 to +3. It occurs due to the change of the Coulomb potential for the 
decreasing [CuC>2-SrO] distance and the corresponding deformation of interfacial 
CuOe octahedra. As a result of the orbital reconstruction, the local electron con- 
figurations of Cu atoms in the antiferromagnetically ordered interface CuC<2 planes 
can be represented by an alterating sequence of {3c/ 5 ^;3t/ 5 ^} local spin states. 
Such an arrangement leads to an enhancement of the local magnetic moment on Cu 
up to 33% which is also consistent with the corresponding Hund rules of electronic 
orbital arrangement. 

Consequently, at the Cu02-SrO interfaces, two basic reconstruction mechanisms 
are involved in order to compensate the interface polarity: (i) electronic charge com- 
pensation which leads to the hole charge predoping of total amount Nh = 1 with 
holes located predominantly in the O 2p orbitals hybridized with Cu 3c/ v 2_ v 2 in the 
CuC<2 planes; (ii) orbital reconstruction associated with changes of the occupancy 
of Cu 3d 3 -2_ r 2 orbitals so that exactly one of these orbitals becomes unoccupied 
by electrons. In contrast to the unrelaxed interfaces where only the first charge- 
predoping mechanism is responsible for the electrostatic neutrality, the relaxation of 
the neutral supercell leads to a combination of two (charge and orbital) mechanisms 
which become equally important in order to achieve the stability of the system. 
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Fig. 7 Projected O 2d orbital density of states in the C11O2 planes at the interface with SrTiOj 
terminated by SrO (LSDA+t/ studies). The zero of energy is at the Fermi level. The top and bottom 
panels correspond to the case of unrelaxed (A = 2.2 A) and optimized (A = 1.83 A) interfacial 
distances. 

The new reconstruction mechanism also involves a redistribution of the partial 
density of states for the O 2p orbitals in the CuC>2 planes as shown in Fig. [7] In the 
unrelaxed heterostructure the Fermi level is located at a distance of about 0.7 eV 
below the top of the O 2p bands implying high hole densities. The interface relax- 
ation leads to a shift of the Fermi level closer to the top of the 2p bands which leads 
to the decrease of the hole charge. It is remarkable that such a modification of the 
structure near the Fermi level is also accompanied by a significant reduction of the 
energy gap between O 2p and Cu 3d orbitals, an effect which is observed in the total 
DOS, displayed in Fig. [8] Such a reduction of the gap is caused by the location of 
the empty Cu 3d x i_ v i band at a distance less then 1 eV above the Fermi level which 
is a clear manifestation of the orbital reconstruction mechanism in the structurally 
relaxed heterostructure. 

It should be noted that up to now, the growth of high quality YBaCuO films on 
the substrates of structurally compatible transition metal oxides remains a challeng- 
ing task due to their roughness caused partially by ionic compensation of interface 
polarities. The existing difficulties with assembling heterostructures as well as with 
probing the interfacial physics directly make the theoretical calculations a powerful 
alternative tool for studies of oxide interfaces. In fact, the importance of the pre- 
sented results is supported by recent resonant X-ray spectroscopic studies of inter- 
faces between YBa2Cu3(>7 and La2/3Ca,/3Mn03 El for which evidence of orbital 
reconstruction with partial occupation of Cu 3d^i_ r i orbitals has been provided. 
Our theoretical studies identify a possible mechanism of such orbital reconstruc- 
tions. In this context, the presented HLRB project contributes new and fascinating 
results in the rapidly developing field of the physics of transition metal oxide sur- 
faces/interfaces. 
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Fig. 8 Total density of states of the superlattice with C11O2 deposited on SrTi03 , terminated by SrO 
(LSDA+f/ studies). The zero of energy is at the Fermi level. The top and bottom panels correspond 
to the case of unrelaxed (A = 2.2 A) and optimized (A = 1.83 A) interfacial distances. 
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